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I investigate non-Gaussian signatures in the context of tachyacoustic cosmology, that is, a nonin- 
flationary model with superluminal speed of sound. I calculate the full non-Gaussian amplitude A, 
its size /nl, and corresponding shapes for a red-tilted spectrum of primordial scalar perturbations. 
Specifically, for cuscuton-like models I show that /nl ~ 0(1), an d the shape of its non-Gaussian 
amplitude peaks for both equilateral and local configurations, the latter being dominant. These 
results, albeit similar, are quantitatively distinct from the corresponding ones obtained by Magueijo 
et. al in the context of superluminal bimetric models. 
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The physics of the very early universe is a rich investigation field, for there is a plethora of potential models that 
solve, at least partially, the well-known problems of the standard cosmological paradigm. Inflationary cosmology [T] is 
surely the most successful model, but alternatives to it have been proposed in the recent past: pre-big bang cosmology 
[5] , ekpyrotic and cyclic models [3] , nonsingular quantum cosmological models [3] , noncanonical models [5J [BJ , among 
others. All of them lead to a near-scale invariant spectrum, so that a very important question emerges: how can we 
falsify such distinct cosmologies? An answer can be provided by the investigation of their non-Gaussian signatures; 
single-field inflation, for example, predicts a nearly-Gaussian CMB anisotropy field, whereas some noncanonical mod- 
i— i els, for example, can predict large deviations from Gaussianity, as in DBI inflation (see, for example, references [7] 
for large-field polynomial models, and [5]). Therefore, the study of non-Gaussian signatures can be a very important 
^2 tool to rule out many models proposed to cope with the problems of the very early universe. 

In particular, noncanonical Lagrangians lead to models with varying speed of sound, a feature which substantially 
enhances or diminishes their non-Gaussian amplitudes, since they are usually related to terms with a c~ 2 dependence 
£■ — [SHU]; hence, for c s < 1, such terms dominate, enhancing non-Gaussianities as in the case of DBI inflation. For 
superluminal models, c s 3> 1, such terms become subdominant, so that the size of the non-Gaussian amplitude becomes 
/nl ~ 1 [IHIIl], which distinguishes this class of superluminal models from inflation (for which /nl ~ 0(0. 01)) 
[141 115j . In |12j . the authors take into account a minimal and a nonminimal bimetric model, that is, models in which 
£N| a disformal transformation between matter and gravity metric is evocated, and has the form g^ v — g^ u — 

where the coupling B is regarded as a constant in the former, whereas in the latter it can run with (f>. They then 
calculate non-Gaussian amplitudes for the nonminimal bimetric model in the superluminal limit, showing that it 
only mildly depends on the tilt n s — 1 of the power spectrum. In [T2] , the authors take a step further and consider 
non-Gaussianity for a wider class of models without slow-roll and exact scale invariance assumptions, including DBI 
and bimetric models, and study the superluminal limit of the latter. However, despite the generality of their analysis, 
there was still room for another class of superluminal models, which we presented in [BJ and dubbed tachyacoustic. 

Tachyacoustic cosmology is a noncanonical, noninflationary model with superluminal speed of sound, in which 
a nearly scale-invariant spectrum is generated by quantum perturbations redshifted outside of a shrinking acoustic 
horizon. Since this model is noninflationary, and the scale factor is a power law of time, a(t) oc t 1 ^, where e = const 
is the flow parameter associated with the ratio —H/H, we can take a radiative equation of state for the model, 
discarding the need of a reheating period to make the scalar field decay into radiation. Also, as shown recently in 
[ho] . DBI and cuscuton-like cosmological solutions also exhibit attractor behavior. Having this model a nearly scale 
invariant spectrum of perturbations, as mentioned, the next step to take is exactly the analysis of its non-Gaussian 
features, which is the subject of this paper. 

As we shall see, the distinguishing feature of the tachyacoustic model with its cuscuton-like Lagrangian 1 with regard 
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1 A cuscuton is a causal field with iniinite speed of sound, originally proposed in 1171 . From the flow hierarchy devised for arbitrary 
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to the other superluminal models explored in [T3J [T3] concerns the contribution of the term £ 3 in the cubic action 
[9] 1 10] to the full non-Gaussian amplitude, giving rise to an extra term with a linear dependence on the flow parameters 
e and s. This dependence is not present either in DBI models (since the coefficient of ( 3 in the third-order action is 
identically null in this case), or in disformal bimetric model, since by projecting its scalar-field action in the Einstein 
frame one obtains a DBI-like action. Therefore, the tachyacoustic model "inherits" an extra contribution from the ( 3 
term, providing a different value for /nl compared to the other superluminal models. Despite small, such difference 
might be observable in the near future, what is of great importance for falsifying superluminal noncanonical models. 

This paper is organized as follows. In Section [II A| I review the basics of tachyacoustic cosmology, and in Section al B| 
I derive some important resul ts concerning the power spectrum. In Section [ill A| I review the third-order action in the 



fluctuations, whereas in 



IIIB 



I calculate the £ contribution to the non-Gaussian amplitudes, which is absent in other 



models. In Sections |III C and |III D] I discuss, respectively, the size and the shapes of the non-Gaussian amplitudes for 
the cuscuton-like tachyacoustic model. In IV I summarize the main results of this paper. 



II. TACHYACOUSTIC COSMOLOGY 



A. Basics 



One of the crucial problems of the very-early cosmology is the presence of a horizon that prevents a causal ex- 
planation for the observed features of the universe. The inflationary paradigm tackles this problem by means of a 
huge accelerated expansion in the very early universe, enforcing physical scales to be deep inside the horizon and, 
in consequence, in causal contact. However, if light could have travelled much faster in the early universe, all scales 
would have been naturally in causal contact, so that the horizon problem would not be an issue at all. This is the 
basic proposal of the varying speed of light theories (VSL) [IS]. Despite VSL theories fail to explain cosmic structures, 
their essential idea remains a very insightful approach to search for alternatives to the inflationary paradigm. 

Along with the speed of light, another velocity parameter comes into play when studying cosmological perturbations: 
the speed of sound parameter c s . In the standard inflationary scenario, c s is a constant and coincides with the speed 
of light. However, in k-essence models, the presence of a noncanonical kinetic term gives rise to the possibility of a 
varying speed of sound. Also, such models possess two horizons, a Hubble horizon and an acoustic horizon, which 
have independent dynamics; since the acoustic horizon depend on the speed of sound, it can be large or even infinite 
at early times for large values of c si solving the cosmological horizon problem without inflation. In this case, as in VSL 
models, superluminal propagation is again evoked, for c s represents the speed of propagation of the perturbations on 
a nontrivial background; then, a relevant question now arises: are these models causally consistent? 

Fortunately, causal problems are usually model-dependent. As pointed out in references [19] and [20], causal 
problems do not affect general k-essence models or bimetric theories; also, cuscuton models are free from them [17) . 
In particular, it can be shown that the noncanonical structure of a general k-essence model induces an "acoustic" 
metric G^" [10]; hence, this two-metric model possess two distinct local causal structures, the first one being the usual 
light cone (connected with the background metric g^ v of the manifold A4), and the second one connected to the speed 
of sound (thereby the "acoustic cone"). It can be shown that if there exists a time coordinate t with respect to the 
background metric (which is everywhere future directed) which plays a role of global time for both spacetimes {AA, g^v) 
and (A4,G~l) , which guarantees the absence of closed causal curves for k-essence models built from homogeneous 
scalar fields on a FRW background |20j . Since the causal structure for small perturbations is determnined by the 
acoustic cone, we can study superluminal propagation safely, for no causal paradoxes arise. 

Tachyacoustic cosmology is a particular solution of a wider class of k-essence models (see Bean et al. [2TJ), which 
we briefly review below. Consider a general Lagrangian of the form C = C[X, <fi), where 2X — g^d^d^cj) is the 
canonical kinetic term (X > according to our choice of the metric signature) . The energy density p and pressure p 
are given by 

p = C(x,<j>), (i) 

p = 2XC X - C, (2) 



k-csscncc models 1211 we reconstructed an action similar to the one associated with the cuscuton plus a canonical kinetic term, which 
accounts for the "cuscuton-like" term we proposed in [6]. 
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whereas the speed of sound is given by 



d = ^=(l + 2X^X\ (3) 



Px \ C 



x 



where the subscript "X" indicates a derivative with respect to the kinetic term. The Friedmann equation can be 
written in terms of the reduced Planck mass M P = l/y/8nG 
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For monotonic field evolution, the field value <p can be used as a "clock" , and all other quantities expressed as 
functions of </>, for example X = X (0), C = C [X (4>) ,(f>], and so on. We consider the homogeneous case, so that 
4> — y/2X . Next, using 

dt ^ d(f> ^ d<f> ' ^ ^ 

we can re-write the Friedmann equation as the Hamilton Jacobi equation 

, o A / 2 

4> = V2X = -tpLffW), (6) 
i-x 

where a prime denotes a derivative with respect to the field cj>. The number of e-folds dN can similarly be re-written 
in terms of dip by: 



/2X 2M P \H'(4>) 



dN = -Hdt = -^=d<p=^(^^)dct>. (7) 
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We introduce a hierarchy of flow parameters, the first three defined by 

(8) 

zU , _ C x 2M P H' (4>) C' x 
HC X C x H (<p) C x 

where we have used equation Note that using definitions (JsJ , we can write down the equation of state for a 
k-essence fluid as 

where we have used equations ([!]) and |2]) . Therefore, the equation of state parameter w can be written in terms of e 
and, conversely, 

e = \{l + w). (12) 

Hence, a matter-dominated (w = 0) tachyacoustic model has e = 3/2, whereas a radiation dominated one (w = 1/3) 
has e = 2. 

In the case of Tachyacoustic Cosmology, we take all flow parameters to be constant, and then solve the flow hierarchy, 
reconstructing the action at the end. In particular, we find that the Hubble parameter and the speed of sound evolve 
as 

H{4,)=H Q [£j S , (13) 
cs = (14) 
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and taking s — —2s, the resulting Lagrangian has a cuscuton-like form [5] 

£(X^) = 2f(cf>)Vx + CX-V(<j ) ), (15) 
where / ((/>) is the analog of the inverse brane tension, V (4>) is the potential and C is a constant given by 

(10) 

To conclude this review, is important to mention that, as pointed out at the beginning of this section and demonstrated 
in [6], tachyacoustic cosmology is free from causal paradoxes. Also, as the acoustic horizon is given by 

D »"^ < 17 > 

we can be shown that its conformal-time dependency goes like 

Dh oc oc e (l-e- S )JV K T (l-e- s )/(l-e)_ (lg) 

Therefore the condition for a shrinking acoustic horizon, 1 — e — s > 0, is no£ identical to accelerated expansion. For 
e > 1 and s < 1 — e, the expansion is non-inflationary, the Hubble horizon is growing in comoving units, and the 
acoustic horizon is shrinking. The initial singularity is at r — 0, and we see immediately that for the tachyacoustic 
solution, the speed of sound in the scalar held is infinite at the initial singularity, and the acoustic horizon is likewise 
infinite in size. Therefore, such a cosmology presents no "horizon problem" in the usual sense, since even a spatially 
infinite spacetime is causally connected on the initial-time boundary. 

There are two other important parameters to be taken into account in our discussion. They appear as the coefficient 
of the term £ 3 when one expands a general noncanonical action up to third-order in the fluctuations, in the following 
combination: 

A = EML-i^+2A, (19) 
where the parameters E and A are given by [21 HO] > 

E = XC X + 2X 2 C X x = (20) 
A = X 2 Cxx + \x^C XX x- (21) 



The parameter A can also be written as [HI [TTJ 
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where fx and the "kinetic part" ex of the flow parameter e = ex + are respectively given by 

f - 63 f _ X dH 

/x ~3eP ex "^M- (23) 



The E term as presented in the right-hand side of equation ( 20 1 holds the same form regardless of the noncanonical 
model involved, but the A term may vary according to the underlying Lagrangian. For DBI-like Lagrangians, for 
instance, A is given exactly by 10J 



Xdbi = (i - O , (24) 



whereas for an arbitrary noncanonical model expression ( 22 ) is generally employed with the assumption that fx is 
constant [UJ. Again, for DBI-like models, from (19), (20 1 and (24) it is straightforward to show that 



Adbi = 0, 



(25) 
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so that the £ 3 contribution to the action vanishes identically. This fact accounts for the absence of the A^ 3 contribution 
to the three-point function amplitude of DBI models. The same happens to disformal bimetric models, as argued in 
the Introduction, for their action reduce to a DBI-like one when projected in the Einstein frame. 
In the cuscuton-like Lagrangian (15 1, we do not have to worry about A, for it is identically null: 



0. 



so that 



A - ?*± ft _ 1 

^CUSC 9 I ^ 



(26) 



(27) 



hence, the coefficient of the £ 3 -term does not vanish, unlike the superluminal models described above, yielding a A^ 3 
contribution to the amplitude. 



Perturbations 



In [B] we have solved the mode equation for an arbitrary noncanonical model [35] in the case of constant flow 
parameters (as defined in equations (8jTo|). In this paper I will take a slightly different path, using the so-called 
"sound-horizon" time dy = c s dr instead of the conformal time r [llj . just for the sake of comparisons with the results 
obtained in [T5] , 

Taking e, s, s constant, using definitions (8p|, and defining 



7 = e + s — 1, 



it follows that 



and 



a-H/) 1 ^, c s ~(- y y/\ H~{-y)-^. 
We next turn to the second-order action for the curvature perturbation £ [22] : 



So 



APp I d 3 xdr z 2 



C ?(VC) 2 



(28) 
(29) 

(30) 
(31) 



where z is defined as z — a\/2e/c s . From now on, a prime ' indicates a derivative with respect to y; then, defining 

q = y/c^z, (32) 



action (31) becomes 



S 2 = ^/dH ir [^-(V,r; 



next, introducing the canonical variable v = MpqC,, from (33) we derive the mode equation 



k - — )v k = 0. 



After a little algebra we can rewrite the last term in the equation above as 

q" if, r 



V 



where 



e + 2s - 3 
2(e + s- 1) 



(33) 



(34) 



(35) 



(36) 



() 



The scale-invariant limit is achieved when q" /q = 2/y 2 , so that v = 3/2 in this case. Since curvature perturbations £ 
are quantum fields, we expand them in terms of creation and annihilation operators, 



C(V, k ) = Uk (y)ctk + u* k (y)ai 



k- 



(37) 



so that the solution for vu(y) associated with the Bunch Davis vacuum is 



(38) 



where H$p are Hankel functions of the first kind. Next, from the modes introduced in (|37|) and solution (|38|) we find 
that 



uk(y) 



1/2 I — tt / l \ 3/2-1/ 

Cs /TT /—rrll),, , ^7 ( ~ k V 



(l + iifcy)e- tt » 



where we have used the following expansion of Hankel functions for small arguments, \ky\ -C 1, 



HP(-ky) = -i 



. 2 v T{v){-ky)-' J 



[1 + iky + O(ky) 



2] -iky 



(39) 



(40) 



together with the assumption that T(v w 3/2) w t/tt/2, which corresponds to nearly scale invariance, as desired. 
From expression (40) we can also compute the derivative of Uk{y) with respect to y, 



<{y) 



2M P yf^Wi \ 2 J 



(41) 



which is plays an important role in the computations of section III B 
From the definition of the C power spectrum 



/i^ 3 iai 2 



(42) 



we find, 



(e + s- I) 2 2 
2{2n) 2 Mpe 



2 o2^-3 Jj2 



(43) 



where the bar refers to quantities evaluated at sound horizon-crossing, y = k . From the computation of the power 
spectrum we find the expression for the scalar spectral index n s , which is related to v and the flow parameters e and 
s by 



n» — 1 = 3 — 2v 



2e + .s 
s + e-1' 



(44) 



III. NON-GAUSSIAN SIGNATURES 



A. Third-order action 



In order to investigate non-Gaussian signatures, we basically compute the three-point function for the curvature 
perturbations, starting from the third-order action in the ADM formalism [51 [TO]: 



S 3 = M P Z / dtd 3 xl-r 3 



E 1 - t + 2A 



^ + ^(e - 3 + 3c^)CC 2 ~ 2a^adQ){d X ) 



2/(0 



SL 



(45) 
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such expression is valid outside of the slow-roll approximation and for any time-independent sound speed. In (45) the 
dots denote derivatives with respect to proper time i, and x 1S defined as 



(46) 



whereas the elements appearing in the last term of the action ( 45 ) are 

1 



/(O = ^c 2 + ^cc 



i((J//2 -mm+o ^himam. + 2a2// 



(ac)(9x)-^ 2 (^(a i c^-x))],(47) 



SL 



dd 2 x 
dt 



Hd 2 X - ed 2 C 



(48) 



where d 2 is the inverse Laplacian. Note that (48) corresponds to the linearized equations of motion, and can be 
absorbed by a field redefinition 



C^C« + /(Cn). 



Since in the limit y — > one has 



H 

J72 



-y) 



v-3/2 



(49) 
(50) 



it is clear from expression (39) that C — > const on superhorizon scales, so that we can drop the second term in (47), 
and we are left with the field redefinition £ — > £ n + "qC^/ {Ac 2 )\ in the present case, for constant flow parameters we 
get i] — 0, so that we can simply skip the last term in the action (45). 



B. Amplitudes 

I next turn to the calculation of the three-point function which leads to the /nl parameter, following the same 
steps as taken in In the interaction picture, the three-point function is given by 



(C(t, k!)C(*, k 2 )C(*, k 3 )> = -i / dt'dat, ki)C(t, k 2 )C(t, k 3 ), H int (t')} 



(51) 



where H- mt = —L iat is the interaction Hamiltonian derived from the corresponding Lagrangian (45), and vacuum 
expectation values are evaluated with respect to the interacting vacuum |f2) (a through discussion on this point can 
be found in [9]); to is some sufficiently early time. Except from the £ 3 term, as discussed in the Introduction, all the 
amplitudes derived have the same form as the amplitudes evaluated in |13| . Then, I shall evaluate here only the A^ 3 
contribution to the non-Gaussian amplitude, which is derived from 



Mr 



Vend „3 

Ay d 3 x 

— oo+ie 



3™ C s v 



c' 3 , 



(52) 



where we have changed the time variable t with the sound-horizon time variable dy = c s dt/a, and used (26). The 
parameter y cn d < indicates the "end" of the deceleration period with infinite speed of sound (in [11] . it refers to 
the end of inflationary or ekpyrotic phase). Next, substituting contribution (37) and |52] ) into ( pjlj ), and using the 
commutation relations [a(k), a^(k')] = (27r) 3 c)( 3 )(k — k'), we find 



<C(ki)C(k 2 )C(k 3 ))/3 = -{2ir) 3 iS < - 3) (k 1 +k 2 + k 3 )u kl (y cnd )u k2 (y cnd )u k3 (y cnd ) 



yc„d 2 

**w 1 1 

—oo-\-ie 



u* kl {y)'u* k2 (y)'u* k3 (y)' + perm. 



(53) 



Also, substituting ( 20 ) , ( 29 ) , ( 39 ) and ( 41 ) into equation ( 53 ) , we find 



(C(k 1 )C(k 2 )C(k 3 )), 3 = -(2 7 r) 3 l ^(k 1 +k 2 + k 3 ; 



2 15-6^ M 4 £ 3 3/2 

y L s,ond 



n end (1 i ] \2-2v i i9/2-3i/ 

(fci/c 2 fc 3 ) |y cnd | 



e 

x — 
7 



oo-\-ie 



/2 



1-3 WV 



2 9/2-3i/ AKy 



perm. + c.c. . 



(54) 



where K = fci + k 2 + £3. The integral appearing in the expression can be simplified as follows: first, from (30 1 we get 
the following expressions, 



1/2 1/2 

Cs C s,end 

where we have defined the exponents 



(-)'• 


JL 


(-) 


V ycnd / 




Vi/ond/ 



- m H cad ( y ^ 



C 5/2 , Wend 
s.cnd 



a = 3 - 2v = 



2e + s 
e + s-r 



2e - s 
e + s-1' 



next, defining 



2»(p) = Im 



2/cnd 



dy 



oo+ie 



V 

Vend 



(-iy) n e iK y 



(55) 



(56) 



(57) 



we can reduce expression ( 54 1 to 



(c(k 1 )c(k 2 )c(k 3 )), 3 = (2^) 3 ^ 3 )(k 1 +k 2 + k 3 )(A :i A ;2 fc3r 



015-61/ A/f 4 ^2 „2 



rr4 

"end I |9— 6w 

I i/cnd I 



-T 2 (/3)-T 2 (a) 



■,cnd 



perm. + c.c. . (58) 



It is important to stress that we took the imaginary part of integral ([57j) for this is the only part relevant for this 
calculation [TT]. For all modes of interest, i^|y C nd| is a small quantity; then, as discussed in [13], the above choice of 
integration contour picks up the appropriate interacting "in-in" vacuum at \y\ — > 00 and takes care of the oscillating 
behavior of the exponential. 

We can further simplify expression (581 by relating the quantities evaluated at y — y cnd to their corresponding 

(59) 



values at sound horizon-crossing, y = K 



H end = H (i^e„d)~ £/7 



c s {Ky end ) sh . 



so that we use the expression for the £ power-spectrum ( 43 1 to absorb the constants H and c s in ( 58 ) ; the final result 
is 



(C(k 1 )C(k 2 )C(k 3 )> e 3 = (2 7 r) 7 <5 3 (k 1 +k 2 +k 3 )P c 2 s ^3^3, 



where 



and 



. _ 3 e + s — 1 ( kik 2 k 3 
~ ~4 c? {^K^ 



n s -l 



lM-c 2 s I t3 (a) 



(60) 



(61) 



l e ,(p) = - (fcxfefca) 2 (K\y CDd \fl 2 (p). (62) 

Using the same methods the other relevant integrals are given by (they coincide with the ones computed in 13J): 

n.-i 



1 / kik 2 k 3 



A, 



4c 2 V 2X 3 
1 / kik 2 k 3 



C(dQ 2 



8c? 1 2if 3 



(e-3)I ci2 ((5) + 3c 2 s I (i2 (a) , 
[(e-2s + l)I m) 2(l3)-c 2 s I (ia02 (a)] 



where 



= -(K\y end \ 



[2b(p) + KX 1 (p)} £ k\k) - X t (p) ~ J2 tffcf \ , 

i<3 i¥=j ) 



I cm *{a) = (K\y cad \r(^2k^j 



2L 2 (p) + Kl_ x (p) +l {p) ^ hkj + fc 1 /c 2 fc 3 Z 1 (p) 



'<3 



(63) 
(64) 

(65) 
(66) 
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Contributions A^q^q x , ■AdQdxd 2 x an( ^ -^(d 2 C)(<9x) 2 ^° ^ ne ^ un amplitude, as found in |13| . are subdominant, for they 
go to zero as c s — > oo, so that we neglect them here. 

Integrals of the form (57) are convergent for p + n > —2 as y enc j —> [TT], and assume the form 



l n (p) = -(K\y cnd \)- p cos P ^r(l+p + n)^-"- 1 ; 



then, in the limit 



oo, the relevant integrals (62), (65) and (66) assume the form 

\2 



(a) 

J C(ac) 2 ( a ) 



K 3 cos— r(3 + a), 



cos "^-r(l + a) 



(2 + a) 
K 



e ^ 2fc i 



2 ( l + o) 
X 2 



E* 2fc 



■ooefrd + a) ^T + ^E^ + t 1 



if 



(67) 

(68) 
(69) 

(70) 



where a is given by (56). Adding up all the relevant contributions to the full amplitude A = A^ 3 + A^ 2 + A^^)2 



given by (61), (63) ana (64) respectively, taking the small tilt (n s — 1 <C 1) and c s — > oo limits, we find that the full 



amplitude for cuscuton-like models is given by 

Ar^— yea - 



kik 2 k 3 \ n3 1 J 3 (kik 2 k 3 ) 



o E ^ + ^ E k * O TO- E k i k 



K^' 1 3 2K 2 



+ (n s - 1) 



9 (M^) 
4 V ; if 3 



i<3 



(71) 



where we have used expressions ( 68]|70 ) and (73). Note that except for the terms containing an explicit dependence 
on the flow parameters, 7 = e + s — 1, equation (28), expression (71) is identical to the one obtained in equation 



(4.15) of |13) : as pointed out in the Introduction, such cont ribu tions are absent in both DBI and nonminimal bimetric 
models, for they come from the £ 3 contribution in action (45). In the next two sections, I discuss the modifications 



induced by such contributions with regard to the full superluminal amplitude found in 1 13] 



C. Size of non-Gaussian amplitudes, /nl 

The parameter that characterizes the three-point amplitude, /nl, is defined as |14j 



(72) 



following [IT], we define /nl at k\ = k 2 = k 3 = K/3, so that 

Ak 1 = k 2 =k 3 



/nl = 30- 



K 3 



For the cuscuton-like models, substituting (71) into (73), we get 

n s -l 

/ 1 \ 

/nl 



5 / -, \ 35 , . 

— (e + s - 1) H h (n s - 1) 

81 V ' 108 y ' 



5 . 25 

— (e + s - 1) 

54 v ; 27 



(73) 



(74) 



where superscript c stands for cuscuton. As for superluminal disformal bimetric theories, the corresponding size of 
non-Gaussianities /nl is given by equation (3.16) in j!3j . 



/nl ~ 0.28 - 0.04e - (1.19 + 0.08e)(n s - 1); 



(75) 
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then, for a red-tilted spectrum n s ~ 0.96, we find, for e = 0.01 and e = 0.3 (these are the e values they have chosen 
to plot Figure 6), 



/ NL ~ 0.327, /nl~ 0.317, 



(76) 



respectively; in the case of matter-dominated (e = 3/2) and radiation-dominated (e = 2) cuscuton-like tachyacoustic 
models one gets 



f: 



NL 



0.266, /, 



NL 



0.232 



(77) 



respectively. Then, cuscuton-like models also produce non-Gaussianities of the same order of magnitude as of the 
other superluminal models, albeit with different magnitudes, as can be seen from the results (76) and ( |77[ ). Such 
differences come from the fact that the £ 3 contribution, present in cuscuton-like models, is not negligible and possesses 
negative sign (for its 7 coefficient, given by equation (28), is negative for the flow parameters e and s chosen), so that 
it reduces not only the size of the non-Gaussianities, but their shapes as well. 



D. Shapes of non-Gaussian amplitudes 

Following the literature [53] , I discuss next the shape of the amplitude at fixed K, by means of the dimensionless 
ratio Ac 3 ^oo/kik2k^. It is convenient to change the variables &2 an d &3 into the dimensionless ones X2 = ^2/^1 and 
x 3 = k^/ki; then, the amplitude to be studied is -4g s ^oo (1,2:2, £3) /X2X3 in the X2 — X3 plane. Without any loss of 
generality the momenta are ordered such that £3 < X2 < 1; the triangle inequality implies x% + X3 > 1. To avoid 
plotting the same configuration twice, the amplitude is set to zero outside the range 1 — x 2 < £32:2 • 

For the sake of comparison, it is instructive to discuss briefly the results obtained in [T^]. First, the full amplitude 
generated is given by their equation (4.15), 



kik 2 k 3 
2K 3 



e^2 k i + K Y1 k itf k i k . 



3 

K ^ -3 2K 2 ^ '" % 



i i^j i<j i^j 

whose shape is depicted in Figure [TJ for the scale-invariant limit, and in Figure [2] for a red-tilted spectrum. 



(78) 




*2 



FIG. 1: Non-Gaussian amplitude 

Acs^oo (l,X2,x 3 ) /x 2 x 3 for n s = 1 without FIG. 2: Non-Gaussian amplitude for the same 

the £ 3 contribution. model, but with n s = 0.96. 

As for the cuscuton-like models, we see in Figure|3]that the extra contribution coming from the (J 3 term modifies the 
shape dependence. In bimetric models, for instance, the amplitude peaks in the equilateral limit k\ = &2 = &3 for both 
scale-invariant and red-tilted spectra; in the latter, besides, the amplitude also peaks in the local limit fc 3 <C k\ = k 2 , 
which dominates over the equilateral mode. Cuscuton-like models share the same features, as can be seen in the 
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Figure [3J also, comparing Figures |]and|] it is clear that the £ 3 -terms contributes to reduce the magnitude of the 
non-Gaussian amplitudes, specially in the equilateral configuration. 

This effect is stronger in the radiative cuscuton-like model, which is expected, since the combination (281 yields 
7 ~ —2.9 for a radiative equation of state with e = 2 and s ~ —3.9, whereas 7 ~ —2.4 for a matter equation of 
state with e = 3/2 and s ~ —2.9. Hence, bimetric and cuscuton-like models, despite their similarities, can be clearly 
distinguished by means of the different results obtained for their sizes and shapes of non-Gaussianity. 



34 • . - . ■ . * 14- 




FIG. 3: Non-Gaussian amplitude ,4c s ^oo (1, 2:2,3:3) /X2X3 for the cuscuton-like model with n s = 0.96 with a radiative equation 
of state (left) and with a matter equation of state (right). 



IV. CONCLUSIONS 



Among the models proposed to solve the puzzles of the very-early universe, those with noncanonical Lagrangians 
possess an interesting property, the possibility of a varying speed of sound. For subluminal models, as in DBI inflation, 
a small speed of sound can produce large non-Gaussian amplitudes as /nl ~ —100, for example; in the superluminal 
case, for instance, as in the nonminimal bimetric model, large non-Gaussianities are suppressed, and the predicted 
size of non-Gaussianities amounts to be around /nl ~ Such model admits a nearly-scale invariant spectrum, 

which induces distortions to the equilateral shape of non-Gaussianities, and hence making it distinguishable from 
other noncanonical subluminal models. 

The cuscuton-like models proposed in [6] also possess a superluminal speed of sound, with the advantage of being 
also noninflationary. This means that, unlike inflationary models, we can choose values greater than one for the flow 
parameter e, thus eliminating the reheating period for a scalar field with a radiative equation of state. Since they 
lead to a nearly-scale invariant spectrum, it is natural to take one step further and investigate their non-Gaussian 
signatures, which wa s the main subject of this work. After deriving the amplitudes for each term in the third-order 
action, equation ( 45 ) , I found that the £ 3 contribution in cuscuton-like models is nonzero, unlike in the nonminimal 
bimetric model, due to its DBI-type action. To check whether this contribution induces any significant deviation 
from the amplitude of bimetric models, I derived an expression for the size of the non-Gaussian amplitudes, given by 
equation (74|, and found that /nl ~ C(l)> which is of the same order as bimetric models. However, a small difference 
in the corresponding values of the /nl parameter was found (equations (76) and (77)), due exactly to the negative 
coefficient 7 emerging from the £ 3 contribution to the three-point amplitude in the cuscuton-like models. Such extra 
terms reduce the size of non-Gaussianities, as well as the height of their shapes in comparison to bimetric models. 
Also, the equilateral configuration becomes subdominant, whereas the "squeezed" triangle mode becomes relevant. 

Then, cuscuton-like models share similar features with superluminal bimetric models, but they have distinct non- 
Gaussian signatures, which is evident from the comparison of the shapes of their corresponding amplitudes. Due to 
the low values predicted for /nl, which are still well below the upcoming constraints to be posed by Planck satellite 
measurements, it is not yet possible to falsify such models observationally but only those predicting huge non-Gaussian 
amplitude sizes. 
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